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Abstract 

By using the Sloan Digital Sky Survey (SDSS) first year type la su- 
pernova (SN la) compilation, we compare two different approaches (tra- 
ditional x an d complete likelihood) to determine parameter constraints 
when the magnitude dispersion is to be estimated as well. We consider 
cosmological constant + Cold Dark Matter (ACDM) and spatially flat, 
constant w Dark Energy + Cold Dark Matter (FuiCDM) cosmological 
models and show that, for current data, there is a small difference in the 
best fit values and ~ 30% difference in confidence contour areas in case 
the MLCS2k2 light-curve fitter is adopted. For the SALT2 light-curve 
fitter the differences are less significant (< f3% difference in areas). In 
both cases the likelihood approach gives more restrictive constraints. We 
argue for the importance of using the complete likelihood instead of the \ 2 
approach when dealing with parameters in the expression for the variance. 



1 Introduction 

By the end of the last century observations of type la supernovae (SNe la), 
used as standard candles, dire ctly established the acceleration of the universe 
(|Riessl . Il998t iPerlmutterl Il999h . an awesome result, possibly only surpassed by 
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the discovery of its very own expansion, around eighty years ago. They are 
still the backbone for the prevailing ACDM concordance model, which is fur- 
thermore corroborated by the combination of other probes, such as, e.g., cosmic 
micro wave background anisotropics baryon acoustic oscillations and galaxy clus 
tering ([Komats ul l201ll:lEisensteinl . l2005[|Percivall20ld: IVikhlininl . 120091: iMantzT 



2010). 



The sample standard deviation in the inferred, uncorrected absolute B mag- 
nitude of typical sets of SNe la is of the order 0.4 mag (even after the exclusion 
of outliers such as th e overluminous 1991T-lik e SNe la and the underluminous 
1991bg-like SNe la) (iPhillipsl . Il993t IVauehanl . Il995t lHamuvl Il996t IContardoL 
2000). This amounts to a fractional uncertainty in luminosity of —30% and 
(under the assumption of negligible uncertainties in all other quantities) in lu- 
minosity distance of 20%. It is thus clear that SNe la are not exactly prima 
facie standard candles at all. However, even such a scatter, as compared to other 
categories of astrophysical sources (other supernova types, gamma-ray bursts, 
etc), is small and, in conjunction with their typical high peak luminosity (10 36 
W~ 4 x lO 9 ^©), justifies the effort to improve their use as cosmological stan- 
dardizable candles. After some phenomenological standardization recipes (cf. 
Section [5]), the scatter in Mb decreases to around 0.15 mag, which amounts to 
a fractional uncertainty in luminosity of —13% and in luminosity distance of 
7%. 

Already i n the landmark p apers from the two original surveying groups 
( Riessl . 1998 : Perlmutterl Il999t ). a legitimate concern was expressed about the 
presence of not properly accounted for systematic effects. At that time, how- 
ever, there was only a small number (the order of tens) of observed SNe la 
in any of the available samples, so the uncertainties were statistically limited; 
the main task was gatheri ng new, larg e r, un i form data s ets. For the recent 
compilations and surv eys (IWood-Vasevi .120071; iKowalskil 120081 iHickenl . 120091 : 
Conlev . 2011 : Kesslerl 20091 Amanullahl . lioiol) . and the more so for the fu- 



ture ones (|Pan- STARRS! : iDEsT ILSSTi ). we are no longer sample-limited; the 
urgency is again on the physics of the SN la phenomenon and all the aspects 
which impact it: nature of progenitor binary sys tem (single-degenerate versus 
double dege nerate, tardy versus promp t events) ( Dildav . 2010t Havden . 201C : 
Maoz , 2010f ) , properties of host g alaxy ^Sullivanl. l2010t iLampeitl . l20ld : Gupta 



20111 ). mec hanisms of explosion ( Hillebrandtl . 120001 ) . extinction/intrinsic color 



variations fNobilil. 120091 lYasudal . l2010h. K-cor rection and tem plate calibration 



( Nugentl . 120021 iHsiaol l2007tl. f lux calibrati on (|Facciolil . lioiiT ). inh omogeneities 
(peculiar velocit ies (jHuil . 120061 iDavis , 2010h . gravitational lensing (Wang, 2005 



Kronborg . 2010. ). etc). We would also like to call attention to the pa rticularly 



clear, i nformative and up-t o-date generic r eview articles on SNe la by iKirshner] 
( 2010T) . Howell ( 20 1 lh . and Goobai ( 2011 ). Simultaneously with this physical 
endeavor from first principles, we should also exercise our best statistical con- 
sistent practices to analyze and mine the data and to test the robustness of our 
inferences. It is to the latter that our paper is devoted. 

In contrast to the traditional \ 2 approach present in most of the literature, 
in this work we discuss a different approach to parameter fitting based on the 
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likelihood. We call attention to the fact that, when we want to estimate both 
the covariance and the mean of a Gaussian process, the ordinary (uncorrected) 
X 2 approach (or any iterative recipe therefrom, for that matter) cannot be 
straightforwardly applied, lest we might lose a nontrivial term in the objective 
function to be extremized. This is necessary when the covariance itself does 
depend on free parameters of the underlying model. Rece ntly, s imilar criticisms 
to the traditional y 2 met hod were p resented by Kim ( Kim . 201 lh (see also 
( Vishwakarmai l2010f 0. In ( Kiml . 12011 ) a likelihood approach to simulated data 
for the distance moduli was used, and the focus was on the determination of 
the intrinsic dispersion, without reference to a particular light-curve fitter. In 
our work, by using real-data, the problem is reconsidered in the context of the 
MLCS2K2 and SALT2 light-curve fitters. In particular, we point out that, in 
principle, with SALT2 additional difficulties could arise due to the astrophysical 
parameters. We compare the results of the traditional \ 2 treatment to the 
likelihood one. We restrict ourselves to SNe la datasets only, without taking 
into account other important sources of information (such as CMB, BAOs or 
clusters) in order to make our point pristine, by avoiding the masking due to 
any other such probes. 

The structure of the paper is as follows. In Section [3J we briefly remind the 
general procedure to go from the raw data to the final estimated parameters. 
In Section [31 we describe the two most widely used "fitter pipelines" (MLCS2k2 
and SALT2), which use the traditional x 2 approach. In Section 01 we critically 
review the aforementioned traditional approach and present the likelihood-based 
one. In Section [H the main numerical results are shown, for both fitter s, in the 
case of the Sloan Digital Sky Survey (SDSS) first year compilation (jKesslei . 
20091 ). Finally, in Section [6j we end up with some discussions and conclusions. 



2 Light-curve fitting 



The "primary" data of any SN la survey are apparent magnitudes (or fluxes) 
in a given set of filters, at a series of epochs (phases), for each supernova. This 
constitutes an array m^y j7 (/i,y, 7 ), where i labels the supernova, Y labels the 
filter (or band) and 7 labels the epoch in the time series. For a given supernova 
i, observed in a given filter Y, the scatterplot of the points (i 7 ,m^y i7 ), where 
tj is the observed time, is what we will call (a sampling of) the raw light curve. 

As mentioned in Section [TJ SNe la are standardizable candles; this means 
there is a phcnomcnological recipe whereby the raw light curves, after being 
subjected to a transformation by a iV s t-parameter function, furnish a new set of 
so-called standardized light curves; this means the dispersion in the new magni- 
tudes is considerably smaller than in the original input set. The aforementioned 
phenomenologica l rec i pe is not unique at all: there are several light-curve fitters 
in the literature (fjhal . 120071: IWand . 1200.4 1200.4- lOuvl 1200.4 120071 : IConlevl . 12008 



Burns! . 2011 ). Here we will exploit some feature s of the tw o most common ones: 
the Multicolor Light-Curve Shape (MLCS 2k2) (I Jhal 120071) one and the Spectral 
Adaptive Light curve Template (SALT2) (|Guvl . l2007l )~ 
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In a nutshell: 



The MLCS2k2 fitting model (jJhal 120071 ) describes the variation among 



SNe la light curves with a single parameter (A). Excess color variations 
relative to the one-parameter model are assumed to be the result of ex- 
tinction by dust in the host galaxy and in the Milky Way. The MLCS2k2 
theoretical magnitude, observed in an arbitrary filter Y, at an epoch 7, is 
given by 



in 



f, ( = M Y , a + PY , n A + q Y , n A 2 + K Y , Y . t +n + + X$™, (1) 



where Y' G {U, B, V, R, 1} is one of the supernova rest-frame filters for 
which the model is defined, A is the MLCS2k2 shape-luminosity param- 
eter that accounts for the correlation between peak luminosity and the 
shape/duration of the light curve. Furthermore, the model for the host- 
galaxy extinction is X Y ° S ^ = C, Y ',~f{a Y i + b Y >/Ry)Ay, where (V',7 := 
X Y ° S ^/ X y ° s q, and ay, by are constants; as usual, Ay is the V band 
extinction, at B band peak (ay = 1, by = 0), and Ry := Ay/E(B — V), 
the ratio of V band extinction to color excess, at B band peak. Finally, 
Xy™ is the Milky Way extinction, K Y / Yil is the if-correction between 
rest-frame and observer- frame filters, and fi is the distance modulus. The 
coefficients My~, p Y >~, and Qy )7 are model vectors that have been eval- 
uated using nearly 100 well observed low-redshift SNe la as a training set. 
7 = labels quantities at the B band peak magnitude epoch. 

Fitting the model to each SN la magnitudes, usually fixing Ry gives 
A, Ay and to, the B-band peak magnitude epoch. 



The SALT2 fitter (jGuvl . 120071) makes use of a two-dimensional surface 
in time and wavelength that describes the temporal evolution of the de- 
redshifted (rest-frame) spectral energy distribution (SED) or specific flux 
(flux per unit wavelength) for SNe la. The original model was trained 
on a set of combined light-curves and 303 spectra, not only from (very) 
nearby but also medium and high redshift SN la. 

In SALT2, the de-redshifted (rest-frame) specific flux at wavelength A and 
phase (time) p (p = at .B-band maximum) is modeled by 

<j>{p, X;x , Xl ,c) = x [M (t, A) + nMi(t, A)] exp[cC(A)]. (2) 

and does depend, through the parameters Xq,Xi, and c on the particular 
type la supernova. M (t, A), M ^{t, A), and C(A) are determined from 
the training process described in ( Guv . 20071 ). Mo(t, A), Mi(t, A) arc the 



zeroth and the first moments of the distribution of training sample SEDs. 
One might consider adding moments of higher order to Eq. ([2]). 

To compare with photometric SNe la data, the (unredshifted) observer- 
frame flux in passband Y is calculated as 

F Y (p(l + z)) = (1 + z) [ d\' [X'd>(p, X')T Y (X'(1 + *))] , (3) 
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where T Y (A) defines the transmission curve of the observer-frame filter Y, 
and z is the redshift. 

As called attention to above, each SN la light curve is fitted separately 
using Eqs. @ and ([3]) to determine the parameters xq, x%, and c with 
corresponding errors. However, the SALT2 light-curve fit does not yield 
an independent distance modulus estimate for each SN la. As we will see 
in the next section, the distance moduli are estimated as part of a global 
parameter fit to an ensemble of SN la light curves in which cosmological 
parameters and global SN la properties are also obtained. 

In the next section we will discuss how to obtain constraints on cosmological 
parameters using MLCS2k2 and SALT2 output quantities as our data. 



3 The traditional % 2 approach 

The prevailing SNe la cosmological analysis is based on the \ 2 function: 

X 2 :=X T '£- 1 X (4) 

where X := (/j, — fi th (z,&)), fi is the set of distance moduli derived from 
the light curve fitting procedure for each SN la event, at redshifts given by z, 
\i th {z, 0) is the theoretical prediction for them, given in terms of a vector of 
parameters and £ is the covariance matrix of the events. 

As discussed in the previous section, each light curve fitter gives a different 
set of output or processed data, which are not related to the distance modulus 
in the same way. In this work, to construct the traditional \ 2 function of this 
section or the proposed likelihood function of the next section, we consider as 
"data" the distance modulus estimation obtained from MLCS2k2 and SALT2 
processing. 



3.1 The x 2 approach from SALT2 output 

The SALT2 light curve fitter gives three quantities, with corresponding errors, 
to be used in the analysis of cosmology: 



:= -2.5 log 



x / dA'M (p = 0,A')T B (A') 



(5) 



to be interpreted as the peak de-redshifted (rest-frame) magnitude in the B 
band, x%, a parameter related to the stretch of the light-curve and c, related 
to the color of the supernova alongside the redshift z of the supernova. The 
distance modulus is modelled, in this context, as a function of (m* B , x±, c) and 
two new parameters, d :— (a.,0), plus the peak absolute magnitude, in B band, 
Mb- Defining the corrected magnitude as 



™ C B ml \S) 



+ ax\ — /3c 



(6) 
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we can write 

ti(S,M B )=m% TI (8)-M B (7) 
Assuming that all SNe la events are independent, one can rewrite Eq. Q 

as 

xl ALT2 (g, S,M(M B , h)) = f) M*, Mb)~ M*; 0, h)? 

^ a t (5) + cr int 

where N is the number of SNe la in the sample, 9 denotes the cosmological 
parameters other than h, with the present value of the Hubble parameter given 
by Hq = 100ft. km • s _1 • Mpc -1 . The theoretical distance modulus is given by 

fith{z;d,h) = 5]og[T> L (z;0)] + im(h). (9) 

with 

/100c/(km/s)\ 
^o(ft):=51og[ - I. (10) 

The dimensionless luminosity distance (in units of the Hubble distance today), 
T>l, for comoving observers in a Robertson- Walker universe, is given by 

V + z ) (vk) Sinh (V*W/'=o Eiho) dz ') - if ^ko > 0, 
V L {z-0)={ (l + ^j^j-fa', ifO fe0 =0, 

t 1 + ^ (v^) (^=^£=0 Eih) dZ ') > if < ' 

(11) 

where O^o is the "curvature density parameter" (whatever the underlying dy- 
namical gravitational theory), such that it is proportional to the three-curvature, 
and 

E(z;8):=H(z;6,h)/H (12) 

is the dimensionless Hubble parameter. As indicated in Eq. ©, the x 2 function 
depends on the parameters M B and h only through their combination 

M(M B ,h):=M B + no(h). (13) 

It may thus be thought of as effectively directly dependent on only the param- 
eters 6, S and M.. 

A floating dispersion term, <7j n t, "which contains potential sample-dependent 
systematic errors th at have not been accounted for and the observed intrinsic 
SNe la dispersion" ( Amanullah . |2010T) . is added in quadrature to the distance 
modulus dispersion, which is given by 

where <r^ 2i is the contribution to the distance modulus dispersion due to red- 
shift uncertainties from peculiar velocities and also from the measurement itself. 



G 



Following (Kesslcr, 2009), we will model this contribution, for simplicity, using 

(15) 



the distance-redshift relation for an empty universe which gives 

1 + Zi 



In 10/ z i (l + z l /2) : 



with 



a z,i = a spec,i + rfec, ( 16 ) 

where a spe c,i is the redshift measurement error, and a pec = 0.0012 is the redshift 

uncertainty due to peculiar velocity. 

As advocated by some groups (|Astierl . [20061) , minimizing Eq. (J8J) gives rise 



to a bias towards increasing values of a and ft. In order to circumvent this 
feature an iterative method is performed, according to their approach. 
In this iterative method, the x 2 presented in Eq. ([5]) is replaced by 

XSALT2 (0, 5, M) - ^ af(ri)+ a L • 

Notice that, in this expression, rj is not a parameter of the X§alt2- in order to 
obtain the best fit values for the parameters, rf is given initial values and the 
optimization is performed on 6,8 and M.. After this step, rj is updated with 
the best fit value of 8 and the optimization is performed again. The process 
continues until a convergence is achieved, which means that rj does not change 
under the required precision. 

During this process <7i n t is not considered as a free parameter to be optimized, 
being determined rather by the following procedure: start with a guess value 
(usually CTi n t = 0.15). Perform the iterative procedure described above. The 
value of (Ti n t is then obtained by fine tuning it so that the reduced x 2 equals 
unity (with all the other parameters fixed on their best fit values) . The iterative 
procedure is repeated once more with this new value and the final best fit values 
are obtained. It is important to note that the value of er int affects both the best 
fit and the confidence levels of the parameters, since it changes the weight given 
to each supernova in the x 2 [cf. Eq. (fI7|) ]. 



3.2 The x 2 approach from MLCS2k2 output 

The MLCS2k2 light curve fitter is also a distance estimator and gives us directly 
a cosmology-independent estimation of the distance modulus. In this context, 
the analogue of Eq. (j8]) is 

XMLCS2k2l t '7 n ) — 2 , _2 T „2 ^ L °> 

. a i i- <7- nt t o^ z i 

where is the distance modulus dispersion as given by MLCS2k2. 

The procedure to obtain <7i n t is similar to the one described in the previous 
Subsection, however, in this case we use only a subsample of nearby SNe la and 
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not the full one, as for the SALT2 analysis. After setting up the value of <7i n t, 
we minimize the XMLCS2k2 using the full SNe la sample to obtain the best fit 
values for and h. 



4 The proposed likelihood approach 

Considering the SNe la light curve fitting parameters as Gaussian distributed 
random variables, we propose to take as starting point the likelihood 

L= - 1 exp(-X T S~ 1 X/2), (19) 

^(2^det£ PV 11 V ' 

which is related to the \ 2 m Eq. (j4|) by 

£ := -2 In L = x 2 + In dot T, + N ln(2n) . (20) 



Eqs. (|T9j) and (|20|) are the single basis upon which our whole statistical proce- 
dure lies. When the full covariance of the problem is known, minimizing x 2 is 
completely equivalent to minimizing C. However, this is not the case for current 
SNe la observations and neglecting the last but one term in Eq. ((20)) would, 
in principle, lead to a biased result. Our proposal is minimizing the following 
functions for each case discussed in the previous section 

N 

C SA LT2(8,5,M,<T int ) = XsALT2(9,8,M,a iat )+YlMvi(t)+<rL) (21) 

i 

and 



N 

£MLCS2k2(0, h, (Tint) = XMLCS2k2( , h, ^int) + ^ MPi + °iat), ( 22 ) 



where we neglected parameter-independent terms. Xg ALT2 (0, 5, cri n t) and 
XMLCS2k2(^' ^' CT int) are given, respectively, by Eqs. © and ([T5]) now considering 
Uint also as a free parameter. With this procedure, we can obtain directly 
unbiased probability distributions functions for all parameters, including CT; nt 
and 8. 



5 Results 

In this section we compare the results obtained from the \ 2 an d the likelihood 
approaches, as describ ed in Sections [3] and Q] using real data from the SDSS 
first year compilation ( Kessler . 20091) . 



In order to perform the comparison, we considered the following cosmological 
models: 

• ACDM, in which we can write the Friedmann equation, in terms of the 
parameters = (f2 m o,f2fco)> as 

E 2 {z; 6) = Q m0 (l + zf + il k0 (l + zf + (1 - Q m0 - Q k0 ). (23) 
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Figure 1: 68% and 95% confidence level contours, in the plane f2 m o x CIao, f° r 
the ACDM model. We marginalized over all other parameters with flat prior. 
The solid (dashed) lines are the results for the likelihood (x 2 ) approach. The 
black circle (gray square) is the best fit value for the likelihood (x 2 )- Left panel: 
Results for the MLCS2k2 version of the SDSS compilation. Right panel: Results 
for the SALT2 version of the SDSS compilation. 

• FwCDM, described by 9 = (Q m Q,w) and 

E 2 (z; 9) = n m0 (l + zf + (1 - n m0 - n kQ )(l + z)^ 1+w \ (24) 



We chose the ACDM and FtuCDM models to directly compare the best-fit 
and the 68% and 95% confidence contours for the parameters 6, <7i nt and 9, for 
both SALT2 and MLCS2k2 data. Th e best fit valu es were obtained with the 
MIGRAD minim ization of the Minuit ( Jamesl 1975 ) implementation in ROOT 



( Antcheva . 2009) and the probability distributions were obtained with Monte 



Carlo Markov Chains (MCMC). We considered as the probability distributions, 
in the context of \ 2 approach, the following functions: 

V SA LT2{0,8,M) = N SALT2 exp[- x 2 SALT2 (9,S,M,a int )/2], (25) 

V M LCS2k2(d,h) = A f MLCS2fc2 exp [~XMLCS2fc2 ( > ^> cr int)/2] , (26) 

where the normalization factors Nsalt2 and NMLCS2k2 are independent of the 
parameters to be estimated. Note that, for the traditional \ 2 method, CT; n t is 
fixed so the probability distribution does not depend on it. 

In Figs. [1] and [2] we show the confidence contours for the parameters 9 for 
ACDM and FiwCDM models, respectively. For these models, we note that the 
differences between the best fit and the area of the contours, for the \ 2 an d 
likelihood approaches, are more significant when the MLCS2k2 fitter is used. In 
fact, for the SALT2 fitter, the differences are not significant (less than 13%) - 
see also Fig. [U and discussion below. If this is a general feature or depends on 
the models or dataset used has to be further investigated. 

In Fig. [3] we show the confidence contours for the SALT2 parameters 8 
for both ACDM and FwCDM models. We can see that there is no significant 
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Figure 2: 68% and 95% confidence level contours, in the plane Q m o x w, for 
the FwCDM model. We marginalized over all other parameters with flat prior. 
The solid (dashed) lines are the results for the likelihood (x 2 ) approach. The 
black circle (gray square) is the best fit value for the likelihood (x 2 )- Left panel: 
Results for the MLCS2k2 version of the SDSS compilation. Right panel: Results 
for the SALT2 version of the SDSS compilation. 



difference in the constraints for a. For fj we find a bias that is, however, small 
compared to the 68% confidence interval for this parameter. 

In Fig. 2] we show the probability distributions for cr int , given by the like- 
lihood approach, for MLCS2k2 and SALT2 data. The traditional % 2 approach 
gives only one value for this parameter without uncertainty and we represent it 
by the dashed vertical line in the figure. We can see that the discrepancy be- 
tween the value obtained from the x 2 approach and the best fit value obtained 
from the likelihood approach is greater for the MLCS2k2 data. The results are 
incompatible at more than 99% confidence interval, which does not happen for 
the SALT2 data. This can possibly be due to the fact that dint is obtained 
using only a nearby sample in the % 2 approach for MLCS2k2 while such dis- 
tinction is not performed in the likelihood approach. This issue deserves further 
investigation and will be the subject of future work. 

We also allow for a possible variation of the parameters a, /3, A4 and (Ti n t with 
redshift, in the context of the SALT2 data. In order to perform such analysis 
the dataset was divided in redshift bins and the cosmological parameters were 
fixed in the best-fit values obtained from the global fit, then releasing a, f3, A4 
and <Ti n t to be determined in each bin. The results are shown in Fig. [5] We 
found evidence of evolution for the parameter /3, in agreement with ([Kessler . 
2009); furthermore we also found evidence of evolution for <Tj nt , which might 



support the use of a variable CT; n t instead of a constant one. 

I n addition to the SDSS compilation, we also used the Union2 compila- 
tion (jAmanullahl . l201dh and the SNLS third year data (jGuyJ, 120111 ). However 
these references do not provide the covariances between the SALT2 parameters 
(m* B ,xi,c). In order to make a fair comparison with the SDSS data, we ran 
again the analyses for all three samples, this time neglecting the cross terms in 
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Figure 3: 68% and 95% confidence level contours in the plane f3 x a. We 
marginalized over all other parameters with flat prior. The solid (dashed) lines 
are the results for the likelihood (% 2 ) approach. The black circle (gray square) 
is the best fit value for the likelihood (x 2 )- Left panel: Results for the ACDM 
model with the SDSS compilation. Right panel: Results for the FwCDM model 
with the SDSS compilation. 



MLCS2k2 5ALT2 




°int "int 

Figure 4: Distributions for rji nt for the ACDM model from the likelihood ap- 
proach. We marginalized over all other parameters with flat prior. The hori- 
zontal lines depict the 68%, 95% and 99% confidence intervals. The dashed line 
is the result for the x 2 approach. Left panel: Results for the MLCS2k2 version 
of the SDSS compilation. Right panel: Results for the SALT2 version of the 
SDSS compilation. 
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Figure 5: Evolution of SALT2 parameters a, /?, M, and <7i n t with rcdshift for 
the SDSS compilation for model ACDM. The cosmological parameters were kept 
fixed in their global best fit values. The black circles and the gray squares are 
the results for the likelihood and the x 2 approaches, respectively. The error bars 
represent the 68% confidence intervals, marginalizing over all other parameters 
with flat priors. 
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Eq. (fH)) . Without taking into account the covariances, the differences in the 
areas for the SDSS compilation are essentially the same. The main difference 
relies on the fact that, for the FiuCDM model, the constraints are tighter for 
the x 2 approach than for the likelihood one. For the Union 2 compilation, the 
difference in the areas ranged from 13% to 14%, while the SNLS data showed 
the greater discrepancies, reaching up to ~53%. 



6 Conclusion 

In this work we considered the SDSS first year compilation and compared two 
different approaches (traditional \ 2 an d likelihood) to determine constraints 
from SN la processed by two of the most used light curve fitters in the literature, 
the MLCS2k2 and the SALT2. The MLCS2k2 gives a cosmology-independent 
estimation of the distance modulus for each SN la, with its corresponding vari- 
ance, which can be directly compared with the model prediction for this quan- 
tity. The SALT2 is not a distance estimator, consequently, we can only have 
an estimate of the distance modulus depending on parameters to be obtained 
simultaneously with the cosmological ones. Furthermore, in the SNe la anal- 
ysis, it is common to introduce a residual, unknown, contribution <7i n t to the 
variance which, in the traditional approach, is determined by imposing that the 
reduced \ 2 be unity, when considering the full sample, in the SALT2 case, or 
only nearby SNe la, in the MLCS2k2 case. 

By comparing the results obtained from the traditional \ 2 approach with 
the likelihood ones, we showed that, for current data and chosen cosmological 
models, there is a small difference in the best fit values and confidence contours 
(~ 30% difference in area) (cf. Figs. □ and © in case the MLCS2k2 fitter 
is adopted. For SALT2 the difference is less significant (< 13% difference in 
areas). In both cases the likelihood approach gives more restrictive constraints. 
We can understand these results by observing, from Fig. 0] that the estimated 
value of (7i n t in the traditional approach, is higher than the peak of the er^t 
distribution obtained with the likelihood method. In fact, for MLCS2k2 the a- m t 
value obtained with the traditional \ 2 approach is outside the 99% confidence 
level of the distribution obtained with the likelihood method. For SALT2 it is 
of the order of 68% confidence level and this might explain why for SALT2 the 
differences between these two approaches are less significant. We should also 
remark that we noticed that the covariance between the SALT2 parameters 
(m* B ,xi,c) has an important role in the above result since we obtained, for 
the 68% confidence level contour, an area ratio of 1.13, when considering the 
covariance, and of 0.95, when neglecting the covariance. 

We also studied the possible evolution of the SALT2 parameters a, j3 and 
Ai with the redshift, splitting the samples in redshift bins and performing the 
fit separately for each one. In this situation, we found evidence of evolution in 
the parameter (3 and also in <7i nt . 

While this paper was in preparation two articles addressing the iss ue of 



using x with unknown variances appeared in the arXiv. The first one (|Kim 
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20111) focused on the question of the determi nation of g; nt , using simulated data 



for the distance modulus. The second one (|Marchl . l201lh studied the SALT2 



case and proposed a sophisticated Bayesian analysis. We did not find in our 
work with the SDSS first year compilation any kind of catastrophic biases of 
parame t ers w hen adopting the likelihood approach, a possibility suggested in 



(jMarch|,L2011) 



In summary, we used current data to compare the traditional x 2 an( i the 
likelihood approaches to determine best fit and confidence regions from SN la. 
We argued that when the variance is not completely known, minimizing the 
traditional \ 2 is n °t formally equivalent to maximizing the likelihood function, 
since the normalization of the likelihood, assumed Gaussian, is also a function 
of parameters to be determined. We conclude suggesting the adoption of the 
likelihood framework instead of the traditional % 2 one, since it is more straight- 
forward, numerically more efficient and self-consistent. 
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